Offline commissioning and performance of a 
multi-reflection time-of-flight mass analyzer with a new 
configuration 


Wen-Xue Huang***, Yong-Sheng Wang®>*¢, Yu-Lin Tian®*, 
Jun-Ying Wang**, Yue Wang**, Zai-Guo Gan®*°¢, Hu-Shan Xu®°4 


*CAS Key Laboratory of High Precision Nuclear Spectroscopy, Institute of Modern 
Physics, Chinese Academy of Sciences, Lanzhou, 730000, China 
’ School of Nuclear Science and Technology, Lanzhou University, Lanzhou, 730000, China 
“Advanced Energy Science and Technology Guangdong 
Laboratory, Huizhou, 516000, Guangdong Prov., China 
School of Nuclear Science and Technology, University of Chinese Academy of 
Sciences, Beijing, 100049, China 


Abstract 


To precisely measure atomic masses and select isotopes of interest, a multi- 
reflection time-of-flight mass spectrometer is being built at SHANS (Spec- 
trometer for Heavy Atom and Nuclear Structure). Its mass analyzer with a 
new configuration different from the existed ones has been commissioned suc- 
cessfully with a chopped !33Cs beam by using a Bradbury-Nielsen gate. The 
potentials on all electrodes have been optimized by using two-dimensional 
potential scan method. With the optimized parameter set, an efficiency of 
50% has been obtained. The mass resolving power with a time-of-flight of 
12.9 ms has been measured to be 74,000, and it has reached 90,000 with a 
longer time-of-flight of ~20 ms. 
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1. Introduction 


In recent years, more than a dozen of multi-reflection time-of-flight mass 
spectrometers (MRTOF-MS) have been or will be constructed all over the 
world [IHI]. The tasks of these devices are to perform accurate mass mea- 
surements and isobaric separation, and fruitful experimental results have 
been achieved already [12H19]. Plaß et al. reviewed the MRTOF-MS de- 
velopments for the research with short-lived nuclei and different instrumental 
implementations. At the Institute of Modern Physics, Chinese Academy of 
Sciences (IMP/CAS), we are developing an MRTOF-MS to attach to the 
end of SHANS (Spectrometer for Heavy Atom and Nuclear Structure [21)) 
to precisely measure atomic masses and select isotopes of interest. 

The conceptual design of the low energy mass spectrometer at SHANS 
was described in Ref. [22]. It consists of a degrader, cryogenic gas catcher, 
quadrupole ion beam guide, radiofrequency quadrupole cooler and buncher, 
MRTOF mass analyzer and time-of-flight detector. The MRTOF mass ana- 
lyzer plays an important role as the ions fly in it and are separated in time 
according to their mass-to-charge ratios. The electrode design and the poten- 
tial applied must be chosen carefully such that the analyzer is time focusing, 
i.e. that time-of-flight differences because of initial energies, positions and 
angles of the injected ions are reduced as small as possible. 


Since the number of electrodes affects the potential smoothness along the 


optical axis and determines the number of power supplies, we designed a new 
type of MRTOF mass analyzer with a different configuration, in which each 
mirror consists of five electrodes rather than four or eight [25]. 
The design and computer simulation have been reported in previous papers 
[26]. 

In this paper, we will concentrate on the offline commissioning of the 


MRTOF mass analyzer and report the achieved performances in detail. 


2. Experimental setup 


Fig. [|shows schematically the experimental beam line used for the offline 
commissioning. A continuous !°3Cs'* beam is produced by a heated alkali 
ion source emitter (a) from HeatWave Labs [27] followed by a skimmer (b) 
and an extraction electrode (c). The emitter is connected to a supporting 
electrode floating at a potential of ~190 V above the ground, thereby defining 
the initial kinetic energy of the ions. The potentials on the skimmer and the 
extraction electrode are 155 V and 100 V, respectively. Then a Bradbury- 
Nielsen gate (BNG, e) [28], which is made of stainless steel wires with a 
diameter of 50 wm and a spacing between wires of 250 um, is employed to 
bunch the beam by switching on/off the +28 V potentials applied on two 
sets of the parallel wires, named as BNG1 and BNG2. 

In order to reduce the emittance by limiting the radial spread of the beam 
and thus improve its quality, three apertures are used in combination. The 
first one is in the skimmer with a diameter of 3 mm, the second one is in the 
flange (d) before the BNG with a diameter of 6 mm and the third one (f) 
is behind the BNG with a diameter of 2 mm. Afterwards, an einzel lens (g) 
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Figure 1: Schematic drawing of the offline commissioning beam line: a— ion source emitter, 
b — skimmer, c — extraction electrode, d & f — collimator, e — BNG, g — einzel lens, h — 


Acc-PDT, i- deflector, j — mirror electrode M1, k — M2, 1— M3, m — M4, n - lens electrode 


L, o — drift tube DT, p — detector. 


with potentials of 0, 35 V, 0 is used to focus the pulsed beam to the middle 
of the acceleration pulsed drift tube (Acc-PDT, h), where the ion bunches 
are accelerated to a kinetic energy of 1500 eV by switching its potential from 
0 to ~1310 V to match the requirement of the MRTOF mass analyzer. A 
cylindrical deflector (i) cut equally into 4 parts is placed between the Acc- 
PDT and the mass analyzer to adjust the ion beam in the radial direction, 
and the potentials applied on these electrodes are 0, 28 V, 38 V and 0, 
respectively. 

The MRTOF mass analyzer consists of four pairs of mirror electrodes 
M1-M4 (j, k, 1, m), a pair of lens electrodes L (n) and a drift tube DT (o). 
The details of the MRTOF structure design can be found in [8] [26]. The 
mirror electrodes, from the outermost ends to the middle, have lengths of 
20, 16, 26, and 26 mm, respectively, the lens electrode has a length of 46 
mm, and the intermediate drift tube has a length of 400 mm. The adjacent 
electrodes are insulated by 4 mm thick ceramic spacers, and the total length 


of this analyzer is 708 mm. The diameter of the electrodes is 60 mm, four 


reduced apertures with diameter of 30 mm located on the outer ends of M1 
and DT are added for electric shielding. 

The ions are captured by using the in-trap potential lift method [29], 
in which the kinetic energy of the ions are reduced below the barrier of 
mirror electrodes by switching off the potential on the drift tube. After 
flying a certain number of laps, the bunches are released by switching on 
the potential on the drift tube and hit a MagneTOF detector from ETP Ion 
Detect™ [30] placed at the end of the beam line. The time of flight (TOF) 
from the BNG to this detector is recorded by a SR430 multichannel scaler 
(MCS) from Stanford Research Systems [31]. 

Two sets of vacuum pump from Edward Vacuum [32] are used to achieve 
an environment with a residual gas pressure better than 3 x 1078 mbar. 
High voltage (HV) modules with high precision and stability from ISEG 
Spezialelektronik GmbH [83] are used to drive the electrodes. Fast solid 
switches from BEHLKE Power Electronics [84] and CGC Instruments [35] 
are employed to switch on/off the voltages varying with time. 

The whole system is controlled remotely by a specialized distributed real- 
time control system based on EPICS [86]. Fig [2{a) shows the fast control 
diagram. A PCI-7813R field programmable gate array (FPGA) module from 
National Instrument followed by a TTL line driver is used to trigger 
the BNG, Acc-PDT, DT switches and the SR430 MCS. The measurement 
cycle is depicted in Fig [2{b). It consists of five channels synchronized with 
a precision of 25 ns. A cycle starts when the BNG is opened. The opening 
time was optimized to be 0.15 us. Twenty three ps later, the potential on 


the Acc-PDT is switched on when the bunched ions reach the middle of the 
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Figure 2: Top: Fast control diagram used for the offline commissioning. Bottom: Overview 


of a measurement cycle (not in scale). 
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Acc-PDT and it lasts 1 ms, which is long enough to ensure that all the ions 
leave the Acc-PDT and be accelerated. The ions are captured 40 ps after 
being bunched when they reach the middle of DT and released after flying 
a certain number of laps, which determines the duration. As the total TOF 
of the ions is much longer than the maximum that the SR430 MCS can be 
recorded, the trigger for the MCS has to be delayed to let the whole TOF 
spectrum be recorded. This delay time is varied according to the number 
of laps which the ions fly in the analyzer, as well as the bin width and the 
bins per record of the MCS. Typically the bin width of 40 ns or 5 ns (the 


narrowest available) are used during the offline commissioning. 


3. Optimization of the MRTOF mass analyzer 


The goal of the optimization process is to find the electrode potentials 
that maximize the resolving power with a good symmetry of the TOF peak. 


The mass resolving power is defined as 
MRP =TOF/2/ATOF, 


where TOF is the total time-of-flight and ATOF is the full width at half 
maximum (FWHM) of the peak. 

Due to the impact of image aberration, field imperfection and mechanical 
misalignment, the peaks of TOF spectrum are asymmetrical [88]. In order to 
quantify the symmetry of the TOF peak and then easily pick out the poten- 
tial parameters with good symmetry during the optimization, the coefficient 


of determination (COD) of Gaussian fit has been calculated. The COD is 


expressed as 
n 


Re =1-— Soy = f(x)? Yi» 


i=1 
where y; is the experimental data and f(x;) is the fitted value with the 
statistical model. The minimum possible value of R? is 0 and the maximum 
is 1. The COD measures how well a statistical model predicts an outcome. 
Since Gaussian distribution is one of the main components of the TOF peak 
and distinctly symmetrical, we can employ the COD of Gaussian fitting to 
describe its symmetry. The closer the COD of Gaussian fitting is to 1, the 
closer the TOF peak has the Gaussian distribution, and the more symmetrical 
the peak is. One cannot determine the range of COD of Gaussian fitting for a 
reasonable symmetry, but after analyzing the spectrum, we found that long 
tails appears in the spectrum when R? is less than 0.980. Note that the 
Gaussian fitting is used if not specified in this paper. 

A Python code was developed to set the timing and potential parameters 
and analyze the experimental results online. The TOF, ATOF and R? 
values were quickly calculated in real time after a measurement finished and 
were stored for future analysis along with the potential parameters. 

The switching times for accelerating ions by the Acc-PDT, capturing and 
releasing ions by the DT were scanned to let these operations start when the 
ions are near the middle of the Acc-PDT and DT. One-dimensional potential 
scans on the skimmer, extraction electrode, BNG and einzel lens and two- 
dimensional on the deflector were performed for narrow beam width and high 
transport efficiency. The optimized values have been described in Section 

For the optimization of MRTOF mass analyzer, potentials on the adjacent 


electrodes were scanned in a two-dimensional way over a given range and a 


linear progression through electrodes as M1 vs. M2, M2 vs. M3, M3 vs. M4, 
M4 vs. Land L vs. DT was followed during the whole optimization process. It 
would be very time-consuming to find the best point(s) in a large space with 
a small step size, so a global search and a local refinement were employed. 
During the global search, the potential ranges were as large as tens of volts 
and the step size was 0.5 V or larger. Some sub-regions with relatively narrow 
distribution and good symmetry were found after the global search. These 
sub-regions were selected for local refinement with a step size of 0.02 V or 
smaller. 

The simulated potentials by using SIMION were chosen as the start- 
ing points of the optimization. All electrodes of the MRTOF mass analyzer 
have been optimized based on 300 laps. In order to have sufficient statistics 
to fit the TOF spectrum, 10,000 bunches were recorded in each measurement. 

As an example, Fig. [3{a)-(c) show two-dimensional TOF, MRP and 
R? contour plots when varying the potentials applied on the M1 and M2. 
Fig. [3{d) shows the MRP contour plots with R? > 0.980. In this case, the 
maximum MRP with a good R? located at M1 = 1073.5 V and M2 = 886.0 
V is the optimal choice. We adopted 1073.5 V as our optimal potential for 
M1 and continued two-dimensional scan on the potentials on M2 and M3. 
This process repeated as needed until all the potentials on the electrodes had 


been optimized. 


4. Results and discussion 


Table |1| lists the optimal potentials applied on the mirrors electrodes, 


lens electrode and drift tube after the whole optimization process. With this 
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Figure 3: (a)-(c): Two-dimensional contour plots of TOF, MRP and R? as a function of 
the M1 and M2 electrode potentials; (d) Same as (b) but with R? > 0.980. 
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Table 1: Optimized potentials and fitted slopes of the time-of-flight deviation with applied 
potential on mirrors electrodes (M1—M4), lens electrode (L) and drift tube (DT). ppm 


stands for parts-in-million. 


Electrode M1 M2 M3 M4 L DT 


Potential (V) 1073.30 884.85 472.67 608.67 —1852.80 672.30 
Slope (ns/ppm) -4.630(7) -2.254(2) 0.670(2) 0.655(2) -0.897(2)  -0.062(3) 


parameter set, we can investigate the TOF spectra at different number of 


laps and study the influence of power supply stability. 


4.1. Performances at different laps 


Fig. [4]shows the measured TOF, ATOF, MRP and R? at different laps 
which the ions fly in the mass analyzer. As the figures indicate, all these 
values are functions of number of laps. The TOF from the BNG to the 
detector increases linearly with the number of laps with a slope of 42.758(1) 
us per lap and the ions can be stored in the analyzer up to 65 ms or even 
longer. 

The M RP increases as increasing number of laps until ~400 and reaches 
a plateau around 90,000. Fig. [5|shows the TOF spectra obtained at 300 laps 
where the optimization process has been taken and 500 laps where the max- 
imum MRP obtained. The MRP reaches 74,000 at 300 laps corresponding 
a TOF of 12.9 ms, and it reaches 90,000 at 500 laps corresponding a TOF of 
21.4 ms. Shown in Fig. aa), the R? is greater than 0.980 in most cases, but 
it becomes smaller and thus the symmetry of the peak becomes worse when 


the number of laps is far from 300, where the optimization was performed. 
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Figure 5: The TOF spectra obtained at 300 laps where the optimization process was 
taken (top) and at 500 laps where the maximum MRP obtained (bottom). The red curves 


are the fitted results with a Gaussian function. 
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The transport efficiency is defined as the ratio of the total number of 
counts after a certain number of laps in the analyzer to the total number 
of counts obtained when shooting through the analyzer. Fig. 6] shows the 
measured efficiency at different laps. The efficiency decreases as increasing 
the number of laps. It clearly shows that about 40% of the ions are lost in 
the first 10 laps and the remaining ions tend to stay in the analyzer for the 
remaining number of laps. The efficiency keeps at ~50% in the 100-500 laps 
and decreases slight afterwards. At 1500 laps, we still have an efficiency more 
than 40%. The losses can be explained as two components, the quick one 
could be due to the acceptance of the analyzer, the emittance of the beam 
or a mechanical misalignment and the slow one due to the collision with 
residual gas inside the analyzer. The relatively higher efficiency compared to 
the values in Ref. should be due to the larger inner diameter and thus a 


larger acceptance of our analyzer. 


4.2. Effect of potential inaccuracy 
Fig. |7|shows how the TOF, ATOF and MRP change when varying the 


potentials on each pair of electrodes. The potential deviation is defined as 
(V —Vopt)/Vopt; Where V is the potential applied on electrode and Vo» is the 
optimized potential listed in Table|1} For all the electrodes of the analyzer, 
the TOF changes linearly with the change of potentials and the sensitivity 
was determined quantitatively by the slopes of linear fit listed in Table |1| It 
shows that TOF is most sensitive to the potentials on M1 and M2, and a 
voltage deviation of 0.1 V on M1 and M2 can result in 0.46 us and 0.25 us 
TOF change, respectively. For MRP, the potentials on M1 and DT are the 


most sensitive ones. It clearly shows that the potential on M1 is most critical 
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Figure 6: The measured transport efficiency as a function of number of laps with the 


optimized parameter set in Table [I] 


to the performance of the mass analyzer. It is understandable that the ions 
turn around at the place of M1 and a small deviation from the ideal potential 
may result in a significant change of the ion trajectory. The change of MRP 
from the DT potential can be explained by the change of the energy of the 
ions trapped inside the analyzer. The MRP is not sensitive to the potentials 


on other electrodes except M1 and DT. 
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Figure 7: The change of TOF (top), ATOF (middle) and MRP (bottom) as a function 
of deviation of potentials applied on the sos ik ppm stands for parts-in-million. Note 


that the points of M3 and M4 are almost overlapped and need to be carefully distinguished. 


4.3. Long time stability test 


In order to observe the behavior of the TOF spectrum as a function of 
elapsed time, we took a long time measurement with the optimized potentials 
listed in Table [I] at 300 laps. Fig.[8]shows the measured TOF and MRP 
obtained in a continuous 52-hour measurement. It clearly shows that the 
TOF drifts with the elapsed time with a period of ~24 hours, which should 
be explained by the drift of the electrical power during the day and night, 
and also maybe the change of the environmental temperature. The maximal 
changes of TOF reached 0.2 us. The MRP also fluctuates, but does not 
have a regular period. 

The symmetry of the TOF distribution can be studied by the 1167 spectra 
obtained during this 52-hour measurement. Fig. |9/shows the sum peak after 
re-centering the TOF peaks to correct for temporal drifts and then adding 
them together. As it shows, the peak is symmetric, except for a small long 
tail on the right side. That could be the result of having a non-ideal trapping 
potential in the analyzer and an asymmetric distribution of the beam at the 
entrance of the analyzer. The sum peak can be fitted very well by using a 
Hyper-EMG function and the ATOF and R? were derived to be 67.9(1) 
ns and 0.999, respectively. The corresponding values by using a Gaussian fit 


are 88.9(1) ns and 0.989. 


5. Summary 


A multi-reflection time-of-flight mass analyzer with a new configuration 
has been constructed and successfully commissioned. The potentials on 


all electrodes have been optimized by using two-dimensional potential scan 
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Figure 8: The measured TOF (top) and MRP (bottom) as a function of elapsed time of 


the measurement with the optimized parameter set at 300 laps. 
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Figure 9: The sum peak of 1167 re-centered peaks obtained in a continuous 52-hour 
measurement at 300 laps in linear scale and in log scale (top right corner). The TOF peak 


is fitted by using a Hyper-EMG function. 


19 


method. 

The coefficient of determination of Gaussian fitting has been calculated 
and used as a probe to quickly and accurately pick out the potential param- 
eters with good symmetry. With the optimized potentials applied on the 
electrodes, the ions can be stored in the analyzer up to 65 ms or even longer. 
An efficiency of 50% has been obtained. The mass resolving power at 300 
laps with a time-of-flight of 12.9 ms has been measured to be 74,000, and it 
can reach 90,000 with a longer time-of-flight of ~20 ms. The mass resolving 
power is limited by the rather broad beam chopped by the Bradbury-Nielsen 
gate and will be improved by using high quality pulsed beam with narrow 
bunch width and small emittance. It will be installed at SHANS to directly 
measure the masses and separate the isobars of fusion evaporation reaction 


products. 
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